Quantum computing of molecular magnet Mni 2 
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Quantum computation in molecular magnets is studied by solving the time-dependent Schrodinger 
equation numerically. Following Leuenberger and Loss (Nature (London) 410, 789(2001)), an ex- 
ternal oscillating magnetic field is applied to populate and manipulate the spin coherent states in 
molecular magnet Mn.12. The conditions to realize parallel recording and reading data bases of 
Grover algorithsm in molecular magnets are discussed in details. It is found that an accurate dura- 
tion time of magnetic pulse as well as the amplitudes are required to design the device of quantum 
computing. 



Quantum phase was proposed to store information in 
connection with a new class of computational algorithms 
based on the rules of quantum mechanics rather than 
classical physics.El One example^s the database search 
problem proposed by Grover .BB Differing from other 
quantum algorithms, the superposition of single-particle 
quantum states is sufficient for Grover's algorithm. The 
Grover's algorithm was successfully implemented using 
Rydberg atoms Jj Recently, Leuenberger and Loss theo- 
retically proposed that molecular magnets Mni2 can be 
used to realize the Grover's algorithm by utilizing multi- 
frequency coherent magnetic radiation.El In the S-matrix 
and time-dependent high-order perturbation theory, they 
showed that it was possible to populate and manipulate 
spin coherent excited states by applying a single pulse of 
weak oscillating transverse magnetic field with an appro- 
priate number of matching frequencies, and the coherent 
state can be applied for storing a multi-bites informa- 
tion. In this paper the population and manipulation of 
a spin coherent excited state in an external oscillating 
magnetic field is studied by solving the time-dependent 
Schrodinger equations numerically. The conditions for 
the field to implement the Grover's algorithm are dis- 
cussed in details. We find that an accurate duration 
time of magnetic pulse as well as the oscillating frequen- 
cies and amplitudes are required to design the device of 
quantum computing in the molecular magnets Mni2- 

Following Leuenberger and Loss, we consider a molec- 
ular magnet with spin S (> 1/2) in presence of a weak 
oscillating transverse magnetic field. The Hamiltonian 
for the system reads 



where 



H — H sp i n + Vi ow (t) + Vhigh(i), 



H spin = ~ASl - BSt + g^ B SH z S z 
V\aw(t) = gn B H Q (t) cos(u> t)S z , 



(1) 



(2) 
(3) 



s-l 



v high (t)= 



gfi B H m {t) 



m—m,Q 



(4) 



H S pi n lis a single spin Hamiltonian for molecular magnet 



Mni2.u | to) are the simultaneous eigenstates of H sp i n and 
the z-component spin operator S z with the energy e m and 
the moment mfi. V\ ow and Vhi g h are the Zeeman terms 
to describe the coupling between the external magnetic 
fields and the spin S. The low frequency is applied along 
the easy axis, and V\ ow supplies the necessary energy for 
the resonance condition. The high frequency transverse 
fields are introduced to induce the transition from the 
initial state \s) to virtual states \m) , to = to , • • • , s — 1, 
and the frequencies uj m , to = too, • • • , s — 2, mismatch the 
level separation by u)q, that is , fru> m = e m — s m +\ + huio, 
and uj s -i — e s — e a +i + ( s ~ m o ~ 1)Tiloq. The phases 
$ m = Y^T=s-i ®k + fm {<fim is the relative phase) and 



H m {t) 



H m , if te (-f, 

0, otherwise. 



) for to = and to > mo, 



The time-dependent Schrodinger equation reads 
8 



ih-\*(t)) = H\*(t)) 



where 



!*(*)> 



E 

771— TTIq 



a m (t)e 



-ie m t/h 



(5) 



(6) 



The eigenstates of to < too have been neglected. The 
molecular magnets Mni2 behaves like single spin, and 
has a ground state with spin s = 10. 

The Grover scheme becomes adapted to describe the 
quantum computational read-in and decoding of the 



1 



quantum data register a = (a s , a s -%, a mo ). Firstly, as 
a simplest example, let us discuss the case of mo = 9. 
At the first stage of the so called "read-in" , the initial 



values of {i 



,(o) 



,(o) 



= 10, 9} are set to be {i 



,(o) _ 



1. 



0}. After irradiating the molecular magnets Mn 12 
with a coherent magnetic pulse with the relative phase 
in the duration time T, one has the quantum data 
register a = (a^, a^) with 

/ gn B H 9 



l io 



cos 



V 2h 



-(9|S_|10)T 



o {1) ■ 

9 sin 



^ 2h 



(9|S_|10)T 



{<) 



(8) 



At the second stage of the so called "decoding" , the so- 
lutions of (^) and (||) become the initial values of {a m } 

in Eq.(||), and the relative phase is set to be ip^ = 0. At 
the end of decoding, one obtains 



,(2) _ 



3 in(2£^(9|S_|10)T) , for 4 



(i) 



0, 



r (!) 

IOr (fig = ir 



(9) 



If the condition of gfi B H g {9\S^\10)T/2h < 1 is sat- 
isfied, the result in Eq. (^|) is in a good consis- 
tence with Ref.[5]. For a longer duration time T ^> 
[gfj,BHg(9\S- 1 10) /27i] 1 ,the condition fails to implement 
the Grover algorithm. 

For mo < 9, it is difficult to solve Eq.([|) analyti- 
cally. However, as it is a single particle problem only 
relating to a finite number of states, it is possible for us 
to solve the differential equations numerically instead of 
evaluating the S-matrix perturbatively as Leuenberger 
and Loss did. The results of explicit numerical calcu- 
lation should be more reliable than other approximate 
results. For illustration, we focus on the case too = 5. 
The parameters used here are the same as the ones in 
Ref.[5]: lu = 5 x 10 7 s~\ T = 10- 7 s, H = H 9 = 2G, 
H 8 /H a = -0.04, H 7 /H = -0.25, H 6 /H = -0.61, 
and H^/Hq = —1.12. For the molecular magnets 
Mni2, A and B in Eq.(|) are 0.56^ and 1.11 x lO^X, 
respectively.Q According to the Ref . [5] , the relative phases 
are ip^ — tp% = ip^ = and ip^ = ip^ = ir for 
encoding the number 13io — H012- Now, we solve six 
first-order differential equations for {a m (t)} in Eq. (|^). 



The initial values of the {a m } are set to be 



and dm = for to ^ 10. 
have aW= (a^j, ag , Or , 



(0) 

'10 



1 



After a duration time T, we 

(1) (1) (1) 



). At the second 

stage for decoding the number, a^ becomes the initial 
values of Eq. (||) and the relative phases should be set by 

(2) (2) (2) n , (2) (2) ^ 

(pg — (p^ — (p 5 = u and Lp s — ip 6 — tt . Unce again, 
the six first-order differential equations with the new ini- 
tial conditions are solved numerically. At the end of de- 



coding, we obtain a( 2 ) = (a 



(2) (2) (2) 
10 ' a 9 ' a 8 



,(2) 



(2) (2)x 



If 



Grover's algorithm is implemented successfully, accord- 
ing to Leuenberger and Loss, one should have 



; (2) 
'10 



1. 



(2) 



,(!) 



2^7, 



.(2) 



0. 



(2) 



2n, 



(2) 



2?7, 



,(2) 



and rj <C 1. Instead, our 



numerical calculation gives 



and 



■,(2) 



,(2) 



0.95, 
0.01, 

; 0.82, 



,(!) 



(2) 



0.31, 
0.02, 

i 0.59, 



-,(2) 



0.06, a 6 2) 



The values of {om } are complex numbers (not real) 
and we only present their absolute values here. There- 
fore, the numerical solution shows that the encoding 
number likes OOOI2, but is not IIOI2 as expected by 
Leuenberger and Loss at the specific values of applied 
field and duration time T . However, {am ,ttm } depend 
on the duration time T, and one may be still able to reach 
to the encoding number IIOI2 by a specific duration time 
T of pulse. 

To clarify the feasibility of the method proposed by 
Leuenberger and Loss, we study a simple, but nontrivial 
case of too = 8 as a representative example. The key ap- 
proximations in Leuenberger and Loss's method are that 
Sm] s = for j < n{= s — too). The transitions from 
|fc) to |to), i.e. S^ k (to < k < s) and all higher order 

amplitudes Sm yS are neglected. In the case of too = 8, 
the first order perturbation are required to be zero: 



0.01, a 



,(2) 



0.04, 
0.01 

1 0.08, 
s 0.01. 



S, 



(i) 
9,10 



and 



5' 



Si 

(i) 
8.9 



(i) 
8,9 " 



0, or at least 



,(2) 



5' 



(2) 
8,10 



S, 



(1) 
9,10 



? (2) 
'9,10 



? (2) 
'8,10 



). In fact, the quantum am- 



plitudes of the transitions by the perturbative S— matrix 
formula are given by 



c(2) _ 2tt / 

°9,10 — — 1 



, , , — ^— ? (9 5-10) — 

V 2h I ujq x 1 1 1 2-k 



g^B 



V 



(10) 



c (2) _ 2tt fgfi B \ 2 H 8 H, 
^8,10 — .1—1 



V 2h ) 



U>0 



■(8|5_|9)(9|S_|10> 



T 
2^ 

(11) 



S, 



(i) 
9,10 



1 gusH^e 



i'i><) 



'8,9 — 



2fl 



1 g[i B U 8 e l 
i 2h 



-(9|S_|10) 



sin(w T/2) 
^/2 



sin (ujqT 2) 
1 — r 1 — - 

wo/2 



(12) 



(13) 



2 



where $9 = (pg and $s = <3>9 + ips ((pg and (ps are the 

(2) (2) 

relative phase). The first two terms {S 9 10 , Sg { } result 
from the second perturbation contributions with energy 



,(2) 



-■9,10 



3 8,10 



conservation. From the condition 
can deduce the field amplitude H$/H = 0.16. However, 
the lower order terms Sg 1 ^ and S^l , in general, are 
not exactly equal to zero even if they do not satisfy the 
energy conservation. They can be zero only if loqT = 2nl 
where I is an integer. Due to these terms Sg~l and S^l , 
the amplitudes of the a%{t) and ag(t) in the state \^(t)) 
are obviously oscillating with varying T. We have shown 
them in Fig. (1). 

Now, let us assume the parameters: T = 1.0 x 10~ 9 s, 
H = 200G, H 9 = 20G and w = 4tt x lOV 1 . It gives 



S, 



(2) 
8.10 



0.11. Following Leuenberger and Loss, we set 



p 9 = and 



,(2) 



corresponds to {ip *• 



0. Then, the recording number OO2 

r, ip^ = tt}, 01 2 to {ip^ = 

0, = tt}, 10 2 to {ip ( 9 1] = tt, ip^ = 0}, and 11 2 

to {p^ = 0, ip^ = 0}. The results by our numerical 
calculation are shown in the Table 1 (where m = 9, 8 
represent the binary digits 2° and 2 1 , respectively). 
Table 1: The numerical results for the case of toq = 8. 



decoded number 


a (1) 
"10 


4 1} 


a (1) 

"8 


(2) 


(2) 

a 9 


(2) 
a 8 


00 2 


0.99 


0.13 


0.10 


1.00 


0.07 


0.01 


01 2 


0.99 


0.13 


0.10 


0.96 


0.26 


0.04 


10 2 


0.98 


0.14 


0.10 


0.98 


0.04 


0.20 


11 2 


0.98 


0.14 


0.10 


0.94 


0.28 


0.19 



Using the above parameters, the proposal by Leuen- 
berger and Loss can be implemented. However, the val- 



ues {\a 



(2), 
10 I' 



,(2)| 



,(2) 



|} are not constants and will be 



changed with the duration time T of pulse. In a cer- 
tain range, 



c(l) 
J 9,10 



or 



c(l) 
^8,9 



tude to 



An) 



will have the close magni- 
which can result in the deformation of 



the recording, or even in the destruction. In Fig.l wc 



plot amplitudes of 
(1) - and - 



,(2) 



and 



,(2) 



with varying T (for 



amplitudes 



,(2)' 



= tt) which corresponds to 012- The 



0.26 and 



,(2) 



0.04 for a duration 



time T = 10 9 s. The value 0.26 can be considered as a 



digit number 1 and 0.04 as 0. The values 



(2) 

1 V ' 



and 



,(2) 



are varied with the pulse duration T. The first maximum 



of 



,(2) 



0.107 at T ~ 1.16 xlO" 



the second max- 



imum 



,(2) 



0.22 at T ~ 1.66 xl0" 9 s which becomes 



,(2) 



0.19 at T - 1.37 



larger than the first minium of 

xl0~ 9 s. Therefore, we conclude that a rather definite 
duration of pulse should be chosen for a stable and clear 
quantum recording. 

One can straightforwardly extend the above discussion 
to the case of m = 5, but it does not change the qual- 



itative conclusion we drew from the case of niQ = 8. 
The more quantum states are used to store information, 
the higher precision are required to design the magnetic 
pulse. It is clear that the lower order perturbation series, 
j < n , of the S-matrix do not vanish for general dura- 
tion T even the relative transitions do not satisfy energy 
conservation. Therefore, to implement the encoding and 
decoding in the molecular magnets by applying an exter- 
nal oscillating magnetic field, one should be very careful 
to choose the appropriate time duration T as well as the 
matching frequencies and amplitudes of the field. 




FIG. 1. Amplitudes 



; (2) 



duration T (for ip { 9 ' = and ip\' = 7r) 



and 

(i) 



as a function of pulse 
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